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Magnetic Resonant Mode in the Low-Energy Spin-Excitation Spectrum 
of Superconducting Rb 2 Fe 4 Se 5 Single Crystals 
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We have studied the low-energy spin-excitation spectrum of the single-crystalline Rb 2 Fe 4 Se 5 supercon- 
ductor (T c = 32 K) by means of inelastic neutron scattering. In the superconducting state, we observe a 
magnetic resonant mode centered at an energy of hco res = 14meV and at the (0.5 0.25 0.5) wave vector 
(unfolded Fe-sublattice notation), which differs from the ones characterizing magnetic resonant modes in 
other iron-based superconductors. Our finding suggests that the 245-iron-selenides are unconventional 
superconductors with a sign-changing order parameter, in which bulk superconductivity coexists with the 
a/5 x a/5 magnetic superstructure. The estimated ratios of ha> ms /k B T c w 5.1 ±0.4 and ha> ms /2A. w 0.7 ±0.1, 
where A is the superconducting gap, indicate moderate pairing strength in this compound, similar to that in 
optimally doped 1111- and 122-pnictides. 

PACS numbers: 74.70.Xa, 74.25.Ha, 78.70.Nx, 74.20.Rp 

Soon after the discovery of arsenic-free iron-selenide su- 
perconductors A 2 Fe 4 Se 5 (A = K, Rb, Cs), also known as 
245-compounds [1], their unprecedented physical proper- 
ties came to light, such as the coexistence of high- T c super- 
conductivity with strong antiferromagnetism [2, 3]. The 
pairing mechanism and the symmetry of the superconduct- 
ing (SC) order parameter in this family of compounds still 
remain among the major open questions. In the majority of 
other Fe-based superconductors, it is widely accepted that 
the strong nesting between the holelike Fermi surface (FS) 
at the Brilliouin zone (BZ) center and electronlike FS at the 
BZ boundary leads to the sign-changing s-wave (s ± -wave) 
pairing symmetry [4]. This scenario has been supported 
by different experimental methods, such as angle-resolved 
photoemission spectroscopy (ARPES) [5], quasi-particle in- 
terference [6], and inelastic neutron scattering (INS) [7, 8]. 
On the other hand, recent theoretical calculations [9] and 
ARPES experiments [10, 11] on the 245-system revealed 
the absence of holelike FS at the BZ center in the electronic 
structure, implying that the nesting between the hole- and 
electronlike FS sheets is no longer present. Hence, several 
theoretical studies proposed alternative pairing instabilities, 
such as d-wave or another type of s ± -wave symmetry with 
sign-changing order parameter between bonding and anti- 
bonding states [12-14]. As a hallmark of sign-changing 
SC order parameter, several authors theoretically predicted 
a resonant mode in the magnetic excitation spectrum be- 
low the SC transition, yet its precise position in momentum 
space still remains controversial [12, 13]. 

A major complication in treating the 245-compounds the- 
oretically arises from the presence of a crystallographic 
superstructure of Fe vacancies that has been consistently re- 
ported both from x-ray and neutron diffraction experiments 
[15]. This a/5 x a/5 superstructure is closely related to the 
static antiferromagnetic (AFM) order persisting up to the 
Neel temperature, T N fa 540 K [16]. Although most of the 



existing band structure calculations have so far neglected 
the superstructure, several others have pointed out that it 
may have a strong influence on the FS shape [17, 18]. How- 
ever, these pronounced FS reconstruction effects have not 
been experimentally confirmed so far [3, 10, 11]. Such an 
uncertainty in the FS geometry and its nesting properties 
makes it hard to predict the exact location of itinerant spin 
fluctuations in reciprocal space. Moreover, band structure 
calculations in the vacancy-ordered magnetic state result in 
insulating solutions for the stoichiometric 245-compound 
[17]. A possible way to reconcile these apparently contra- 
dictory observations is to assume a nanoscale phase sepa- 
ration of (i) insulating vacancy-ordered magnetic domains 
and (ii) metallic non-magnetic phase domains with effec- 
tive electron doping that could host superconductivity at 
low temperature. Such kind of electronic phase segrega- 
tion, resembling the situation in hole-doped 122-pnictides 
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Fig. 1. The dc magnetic susceptibility measurements on three 
representative single-crystalline Rb 2 Fe 4 Se 5 samples from the same 
batch. A sharp diamagnetic response is observed in the ZFC mea- 
surement right below 32 K, indicating ~100% exclusion of the 
external magnetic field. 
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Fig. 2. (a) The in-plane projection of twinned magnetic and nu- 
clear Bragg peak positions arising from the -/5 x -/5 Fe-vacancy 
superstructure. The two sets of dots and the corresponding dashed 
lines (red and blue) denote magnetic superstructure Bragg re- 
flections and the magnetic BZ boundaries for the right and left 
twin domains, respectively. Black dotted lines represent the Fej 
unfolded BZ boundary. The arrow shows the trajectory of elastic 
momentum scans, (b) Elastic scans along the K direction [arrow 
in panel (a) and equivalent scans] measured at 1.5 K. The almost 
identical neutron intensities of two symmetric magnetic Bragg 
peaks indicate the nearly equal population of the two twin do- 
mains in the sample, (c) INS intensity at 11.5 meV in the vicinity 
of the magnetic ordering wave vector (1.3 0.1 0.5) F . The intense 
spin-wave excitation peak is consistent with recent time-of-flight 
INS measurements on an insulating Rb 2+5 Fe 4 Se 5 compound [26]. 

[19], found recent support from ARPES [3], scanning nano- 
focus single-crystal x-ray diffraction imaging [20] , scanning- 
tunneling microscopy [21], and optical spectroscopy [22]. 
In this letter, we provide experimental insight by using INS 
to directly probe the elementary magnetic excitations in 
superconducting Rb 2 Fe 4 Se 5 (RFS). 

For the present study, we used a mosaic of RFS single 
crystals with a total mass of ~ 1 g, grown by the Bridg- 
man method [23]. The nearly stoichiometric and homo- 
geneous composition with Rb:Fe:Se = 0.796:1.596:2.000 
(1.99:3.99:5) has been determined by wave-length disper- 
sive x-ray electron-probe microanalysis using a Camebax 
SX50 analyzer with an accuracy of 0.5% for Fe and up to 
1% for Se. The SC properties of the sample were char- 
acterized by magnetometry, where ~100% flux exclusion 
was observed in the zero-field-cooled (ZFC) measurement 
for temperatures up to T c = 32 K [Fig. 1(a)]. The INS ex- 
periment was performed at the thermal-neutron triple-axis 
spectrometer IN8 (ILL, Grenoble), with the sample mosaic 
mounted in the (HH0)/(00L) or (H00)/(00L) scattering 



planes. The wave vectors Q = (0.5 0.5 L) Fe and (0.5 L) Fe 
were directly accessible in our scattering planes, and the 
spectrometer further allowed us to tilt the sample in order 
to access Q-vectors in a certain range out of the scatter- 
ing planes. Here and henceforth, we are using unfolded 
reciprocal-space notation corresponding to the Fe sublattice, 
which we denote as Fe l5 because of its simplicity and the nat- 
ural correspondence to the symmetry of the observed signal 
[24]. We quote the wave vector Q = (HKL) in reciprocal- 
lattice units (r.l.u.), i.e. in units of the reciprocal-lattice 
vectors a*, b*, c* of the Fe sublattice (a* = 2n/a, etc.). 
Here a = b = 2.76 A is the room-temperature distance be- 
tween the nearest-neighbor Fe atoms, and c = 7.25 A is the 
distance between Fe layers. All INS measurements were 
done in the fixed-fc f (fc f = 2.662 A -1 ) mode, using double- 
focused PG(002) monochromator and analyzer. A 5 cm 
thick oriented PG filter was installed before the analyzer 
to eliminate higher-order neutron contamination, and no 
collimation was applied to maximize the intensity. 

We start with the magnetic Bragg peak patterns arising 
from the </5 x </S Fe-vacancy superstructure. Panel (a) in 
Fig. 2 is a sketch of magnetic and nuclear superstructure 
Bragg reflections, projected on to the two-dimensional Qm = 
(H K) plane. Black dots and the large dashed rectangle 
correspond to the center and boundaries of the unfolded Fe : 
BZ, respectively. The solid dots represent magnetic Bragg 
reflections from two twin domains, and the corresponding 
dashed lines, rotated clockwise and counterclockwise with 
respect to the ¥e 1 BZ, are magnetic zone boundaries of 
the two twin domains. Forbidden magnetic Bragg peaks, 
which coincide with the nuclear Fe-vacancy superstructure 
reflections seen by x-ray diffraction [15], are shown by 
empty circles. Figure 2 (b) shows elastic scans crossing two 
magnetic Bragg peaks at (0.7 ±0.1 0.5) Fe , as shown by the 
arrow in panel (a), and along equivalent reciprocal-space 
directions in higher BZs. Note that the two magnetic Bragg 
peaks at if = ±0.1 originate from different twin domains, so 
that their similar intensity indicates almost equal population 
of both twins in our sample. The magnetic Bragg peak 
intensity decreases more rapidly when moving to a higher 
BZ along the out-of-plane direction than in-plane, indicating 
that the magnetic moment is oriented predominantly along 
the L -direction in this system. This is consistent with the 
reported spin configuration in the magnetically ordered 
phase, in which spins are alternatively pointing up and down 
along the c-axis [25]. Figure 2 (c) shows inelastic magnetic 
intensity in the vicinity of the AFM ordering wave vector Q = 
(1.3 0.1 0.5) Fe at 11.5 mey measured at low temperature, 
T = 1.5 K. The intense spin-wave peak is consistent with 
recent time-of-flight INS measurements on an insulating 
Rb 2+ 5Fe 4 Se 5 compound [26]. 

Now we turn to the INS measurements across T c near 
a few candidate Q-vectors, where the magnetic resonant 
mode could be expected. Figure 3, (a) and (b), displays 
raw energy-scan spectra recorded above and below T c at 
Q = (0.5 0.3125 0.5) Fe , where the resonance has been 
theoretically predicted [12], and at (0.5 0.5) Fe , where it 



800 -, 

CD 

1 700- 
E 

| 600- 

O 

o 

5 500- 

■(/J 

tz 

CD 

."S 400- 

c 

o 

1 300 

~ 60 

CD 

1 40 
E 

2 20 

I 

o 

S- -20 

IT) 

en 

T -40 
f. -60 



■ (a) (0.5 0.3125 0.5) • 1.5 K S 


(b) (0.5 0.5) • 1.5 K 


- 400 ° as K 


° 35 K 


. 350- 
300- 


,8 

2a»! #, o 

•»o o 


e 


• 

S 

- -Vw*' - 


c 


12 14 16 • 

cP 




v" 




°8SS°i*^ 








" (c) • (0.5 0.3125 0.5) 


"(d) ■ (0.5 0.5) 


> (0.5 0.25 0.5) t Li- T T 


• (0.5 0) 


■ (0.5 0.5 0.5) ijffltij 




ittr- §T\fl 


ij I 


1 Al*Wr +Jw4— 


i .^S^sS^- 


{mi#^ff^ 


pflf!P»^ ? *f 




/r £ 




;/r : 



2 4 6 8 10 12 14 16 18 2 4 6 8 10 12 14 16 18 20 
Energy (meV) Energy (meV) 




320- 

10 20 30 40 50 
Temperature (K) 



0.0 0.1 0.2 0.3 0.4 0.5 
K(r.l.u.) 



Fig. 3. (a, b) Raw energy scans measured in the SC (1.5 K) and 
normal (35 K) states at Q = (0.5 0.3125 0.5) Fei and (0.5 0.5) Fei , 
respectively. The inset in panel (a) shows the zoomed-in part of the 
resonant peak in the raw data, (c) Intensity difference between the 
SC state and the normal state at three Q-vectors: (0.5 0.25 0.5) Fei , 
(0.5 0.3125 0.5) Fei , and (0.5 0.5 0.5) Fei . While there is no pos- 
itive intensity at (0.5 0.5 0.5), a clear resonance peak (shaded 
region) is observed around 14 meV both at (0.5 0.25 0.5) Fei and 
(0.5 0.3125 0.5) Fei . (d) The same plot as in panel (c), but for 
Q = (0.5 0.5) Fei and (0.5 0) Fei , where the magnetic resonant 
mode has been found in other Fe-based superconductors, but is 
absent here. The base line in (c) and (d) is the difference of the 
Bose factors, (e) Temperature dependence of the raw INS inten- 
sity at 14 meV and Q = (0.5 0.3125 0.5) Fei that demonstrates an 
order-parameter-like behavior with an onset at T c . (f) Intensity 
difference of momentum scans along the BZ boundary, measured 
below and above T c , with a maximum at the commensurate wave 
vector Q res = (0.5 0.25 0.5) Fei . The solid line is a Gaussian fit with 
a linear background. Different symbols represent identical mo- 
mentum scans measured in different experiments, rescaled to the 
(002) nuclear Bragg peak intensity. The position of the resonant 
mode predicted by Maier et al. [12] is shown by the arrow. 

is usually found in other Fe-based superconductors [7, 8]. 
In the absence of any resonant enhancement, the intensity 
is expected to be higher in the normal state due to the 
influence of the Bose factor at low energies. Already in 
the raw data, one can see that this is the case for all data 
points except a narrow energy region around 14meV at 
Q = (0.5 0.3125 0.5) Fei . 

To emphasize this effect and to eliminate the energy- 



dependent background, we plot temperature differences of 
the same datasets in Figs. 3 (c) and (d). Also shown are the 
difference spectra for Q = (0.5 0.5 0.5) Fei , (0.5 0.25 0.5) Fei , 
and (0.5 0) Fe . As seen in Fig. 3 (c), a prominent peak 
(shaded region) is found at Hco ies ss 14meV for Q = 
(0.5 0.3125 0.5) Fei and Q = (0.5 0.25 0.5) Fei , which we 
attribute to the magnetic resonant mode. However, no such 
peak is observed at Q = (0.5 0.5 0.5) Fe , in contrast to some 
alternative predictions based on the d-wave pairing symme- 
try [13]. Figure 3 (d) also demonstrates the absence of any 
resonant mode at Q = (0.5 0.5) Fe and (0.5 0) Fe , where 
it is usually found in iron pnictides [7, 8]. At these wave 
vectors, the data simply follow the solid line, which is the 
Bose-factor difference between 1.5 K and 35 K. 

To verify whether the observed redistribution of spectral 
weight at low temperatures is related to the SC transition, 
we have measured the temperature dependence of the reso- 
nance intensity at Q = (0.5 0.3125 0.5) Fe , which is shown 
in Fig. 3 (e). Indeed, an order-parameter-like increase of in- 
tensity below T c is found, which is accepted as the hallmark 
of the magnetic resonant mode. 

To pin down the exact location of the resonance in Q- 
space, we have measured momentum scans along the BZ 
boundary at both temperatures. Their difference is pre- 
sented in Fig. 3 (f) and suggests a maximum at the commen- 
surate nesting wave vector Q res = (0.5 0.25 0.5) Fe shown 
by the star symbols in Fig. 2 (a), close to the predicted res- 
onance position, Q = (0.5 0.3125 0.5) Fei [12]. Yet, the 
disagreement is small compared to the Q width of the peak, 
which explains the similar INS response at both Q-vectors, 
as seen from Fig. 3 (c) . Because the position of the nesting 
vector is strongly doping dependent, and the calculations 
in Ref. 12 were done for the arbitrary doping of 0.1 elec- 
trons per Fe, a quantitative agreement with our results is 
not expected. The observed wave vector suggests an even 
higher level of the effective electron doping in the metallic 
phase of the sample, which is difficult to reconcile with 
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Fig. 4. Normalized resonance energy, Ha> les /k B T c , in Fe-based 
superconductors for q z = and q z = n [30]. This ratio for RFS is 
slightly higher than for 122-compounds, but comparable to 11-, 
111-, and 1111-type superconductors. 



the stoichiometric chemical composition unless we assume 
electronic phase segregation into electron-rich and electron- 
poor regions of the kind discussed in Ref. 20 and 21. Such 
a high doping level of the metallic regions would also agree 
qualitatively with the ARPES results [3, 10, 11]. 

By comparing the normalized resonance intensity in RFS 
with that in the nearly optimally doped BaCFej.^Ni^Asj 
(BFNA), measured in a similar experimental configuration 
at the same spectrometer [24], we find that the intensity 
at the resonance energy in RFS is approximately a factor 
of three smaller than in BFNA. Because we expect four 
nonequivalent resonant peaks in the BZ of RFS from sym- 
metry considerations, as opposed to only two such peaks 
[Qres = (°- 5 0)^! and (° 0.5) Fe J in the 122-system, the total 
resonant spectral weight in both compounds turns out to be 
comparable. 

It has been shown that the resonance energy scales lin- 
early with T c in iron-based superconductors, with a ratio 
of Hco res /k B T c that slightly varies between different families 
[8, 24, 27], but is generally lower than the respective ratio 
of ~ 5.3 measured in high-T c cuprates [28]. In Fig. 4, we 
compare this ratio in all Fe-based superconductors, in which 
the resonant mode has been found [see Ref. 30 and refer- 
ences therein]. The value for RFS, fia> res /k B T c w 5.1 ± 0.4, 
lies slightly above the nearly universal ratio of 4.3 estimated 
for 122-compounds (solid line) [24], but is close to that in 
FeTe 1 _ JC Se x , LiFeAs, La-1111, and cuprate superconductors. 

Another important dimensionless parameter that allows 
an assessment of the pairing strength in unconventional su- 
perconductors is the Ra> res /2A ratio, where A is the super- 
conducting gap. We have recently shown that in Fe-based su- 
perconductors this ratio varies between the universal value 
of ~ 0.64, typical for strong-coupling superconductors with 
high T c [29], and unity for low-T c compounds [30]. In 
the 245-systems, the SC gap has been measured by ARPES 
and NMR [11, 31], producing the average 2A/k B T c ratio 
of ~ 7.2 ± 0.4 [30]. It corresponds to the Hco ies /2A ratio 
of ~ 0.7 ± 0.1 in RFS, slightly above the strong-coupling 
limit and supporting the general trend found in all Fe-based 
superconductors [30]. 

To conclude, we have observed the magnetic resonant 
mode coexisting with the AFM order in RFS below the SC 
transition. Our finding suggests unconventional pairing 
with a sign-changing order parameter in the 245-systems, 
qualitatively consistent with theoretical predictions, made 
under the assumption of finite electron doping in the metal- 
lic phase volume. Although this result points toward an 
electronic phase separation scenario [20], the resonant 
spectral weight is comparable to that in optimally doped 
122-compounds. It evidences bulk superconductivity in our 
sample, consistent with the ~100% flux exclusion observed 
in the ZFC magnetization measurement, and implies an 
unusual coexistence or microscopic segregation of strong 
antiferromagnetism with superconductivity, which should 
be addressed by future theoretical work. The estimated ra- 
tios of Hco res /k B T c k 5.1 ± 0.4 and fiw res /2A w 0.7 ± 0.1 in 
this compound indicate moderately strong pairing, similar 



to other Fe-based superconductors. 
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